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Abstract: 

In this paper, I investigate a local effect of polarization of the Cosmic Microwave Back¬ 
ground (CMB) in clusters of galaxies, induced by the Thomson scattering of an anisotropic 
radiation. A local anisotropy of the CMB is produced by some scattering and gravitational 
effects, as, for instance, the Sunyaev Zel’dovich effect, the Doppler shift due to the clus¬ 
ter motion and the gravitational lensing. The resulting anisotropy A I/I depends on the 
physical properties of the clusters, in particular their emissivity in the X band, their size, 
their gravitational potential and the peculiar conditions characterizing the gas they con¬ 
tain. By solving the Boltzmann radiative transfer equation in presence of such anisotropies 
I calculate the average polarization at the centre of some clusters, namely A2218, A576 
and A2163, whose properties are quite well known. I prove that the gravitational effects 
due to the contraction or to the expansion have some importance, particularly for high 
density structures; moreover, the peculiar motion of the cluster, considered as a gravita¬ 
tional lens, influences the propagation of the CMB photons by introducing a particular 
angular dependence in the gravitational anisotropy and in the scattering integrals. Thus, 
the gravitational and the scattering effects overally produce an appreciable local average 
polarization of the CMB, may be observable through a careful polarization measurements 
towards the centres of the galaxy clusters. 
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1. INTRODUCTION 


The first detection of the CMB anisotropy, performed by COBE (Smoot et al. 1992) 
four years ago, opened a large debate about the presence of a detectable polarization in 
the background radiation: in fact, one expects (Rees 1968) that the Thomson scattering 
of such an anisotropic radiation on the electrons plasma produces a linear polarization. 

The CMB propagation is generally studied by using the Boltzmann transfer equation 
(Peebles 1980); the Chandrasekhar formalism (Chandrasekhar 1960) is particularly useful 
to study a polarization problem as the one I want to handle. After the introduction of 
a scattering matrix that depends on the polar angles (O ', (/>') of the incoming beam, the 
evolution of the radiation with the time, the photon energy and with the space coordinates 
is described by using a partial derivative differential equation, containing, as variables, the 
four Stokes parameters fixing its polarization status at all the times. 

One should distinguish between a global and a local effect of CMB anisotropy and 
polarization: the global anisotropy is produced by some primordial density fluctuations 
(see, for instance, Harari & Zaldarriaga 1993) and by some time-dependent fluctuations 
of the space-time metric (i.e. by gravitational waves, originated, for instance, during the 
inflationary era (Sachs & Wolfe 1967). 

I discussed these effects, particularly those due to tensor perturbations, in (Gibilisco 
1995; Gibilisco 1996); similar approaches to this problem can be found in (Basko & 
Polnarev 1980; Polnarev 1985; Crittenden et al. 1993a and 1993b). 

Here, I want to study the CMB transfer equation from a local point of view, i.e. 
considering the anisotropy effects in clusters of galaxies: in this case, the clusters proper¬ 
ties determine in a fundamental way the characteristics of the CMB photons which cross 
them. In particular, the possible sources of CMB anisotropy in cluster of galaxies are the 
gravitational bound (GB) effect, the collapse (C) or expansion (E) effect, the gravitational 
lensing (GL), the Sunyaev Zel’dovich (SZ) effect and finally a Doppler-like shift (D): all 
these effects will be discussed below. 

In this work I prove the importance of the gravitational phenomena: in particular, the 
motion of the cluster, considered as a lens, introduces an additional angular dependence 
in the scattering integral, thus conditioning the cancellation of some terms there; the 
gravitational collapse effect dominates for massive clusters, having a high density, while 
the Doppler shift is always negligible if compared to the SZ effect; note however that it 
may have some importance in a different framework, namely in galactic halos: in this case, 
one expects a total irrelevance of the SZ and of the gravitational effects, due to the low 
values of density and temperature characterizing the halos; only a significant fraction of 
MACHOS within them could enhance the relevance of the gravitational phenomena, but 
more theoretical investigations are necessary to clarify this problem. 

Here, through the mathematical solution of the Boltzmann transfer equation, I cal¬ 
culate the overall CMB polarization induced by the Thomson scattering of this locally 
anisotropic radiation and I evaluate the relative importance of these effects. In the follow¬ 
ing, for simplicity, I will assume to have a spherical shape for the clusters and an isothermal 
model for the gas they contain (see, for instance, Boynton et al. 1982). 
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This paper is structured as follows: in sect. 2 I will briefly discuss the theory of the 
radiative transfer, by recalling the main equations governing the photons propagation in 
presence of Thomson scattering; I will also recall the formalism of the Stokes parameter 
describing the polarization status of the CMB. In sect. 3 I will discuss in some detail the 
peculiar effects contributing to create a local anisotropy for the CMB crossing the clusters; 
their dependence on the polar coordinates will be examined and the consequences of such 
a behaviour will be pointed out. In sect. 4 I will solve the transfer equation by studying 
the space evolution of the Stokes parameters on time-constant hypersurfaces; then, I will 
show my results for the average polarization degree in some clusters, namely A2218, A576 
and A2163. Finally, in sect. 5 I will present my conclusions, emphasizing the importance 
of the gravitational effects in producing a local polarization of the CMB. 

2. THE THEORY OF THE RADIATIVE TRANSFER 


A detailed discussion of the radiative transfer theory in presence of Thomson scattering 
was given by Chandrasekhar (Chandrasekhar 1960); the fundamental equation governing 
the evolution of the photon distribution function is the Boltzmann equation (for a detailed 
explanation see, for instance, (Peebles 1980)) mathematically expressing the Liouville the¬ 
orem of the conservation of the probability density in the phase space. 

A suitable way to write the Boltzmann equation adopts as variables the four Stokes 
parameters, locally describing the polarization status of the radiation (Chandrasekhar 
1960); these parameters are used as components of a 4-vector n a — n a (9, </>, v), which is a 
function of polar angles ( 0 , <f>) and of the photon frequency v. 


n a = (h, I r , U, V). 


( 2 . 1 ) 


In eq. (2.1) I), I r are the left and right intensities, defined as in the Chandrasekhar 
geometrical formalism (Chandrasekhar 1960) based on the polarization ellipse (see fig. 1); 
with this definition, the total intensity of the radiation is I = Ii+I r , while Q — Ii—I r and 
U are the parameters indicating the presence of a linear polarization; finally, V represents 
the circular polarization, but its evolution can be neglected because it totally decouples 
from the others parameters (see Appendix). As a result, the Thomson scattering of an 
anisotropic radiation produces linear polarization only, as proved by Rees in (Rees 1968). 

The final expression for the radiative transfer is a partial differential equation that 
reads: 

/ dn (y dn \ dv dn ot N e R(r /) 

V dr] ^ dx a ) di] du Atx 

r- rl r27r 


X 


47m + 


n(//, (/)') P([j,, </>, //, ()>') djj! dcj) 1 


( 2 . 2 ) 


F -i 


In eq. (2.2), r/ is the comoving time, defined as f[dt/R(t)\, and R(t) is the scale factor of 
the Universe; 7 0 are the components of an unit vector in the propagation direction of the 
photons, or = 6.65 x 10 -25 cm 2 is the Thomson scattering cross-section, /r = cos 6 and 
P ( /_i. 0, [j!. cfR) is the Chandrasekhar scattering matrix, whose explicit form as a function of 
the polar angles is given in (Chandrasekhar 1960) and partially shown in the Appendix; 
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finally, N e is the comoving number density of free electrons, depending on the ionization 
history of the Universe. 

An analytical way to solve the transfer equation in presence of cosmological gravi¬ 
tational waves has been discussed in (Gibilisco, 1995) on the basis of the theory of the 
Volterra integral equations: here it is sufficient to recall that eq. (2.2) contains a col¬ 
lisionless part, given by the time component of the geodesic equation, and a collisional 
term, represented by the integral in its right-hand side. That really gives all the informa¬ 
tions about the scattering processes, responsible for the production of a linearly polarized, 
outgoing radiation. 

3. LOCAL EFFECTS PRODUCING THE CMB ANISOTROPY AND 
POLARIZATION IN CLUSTERS OF GALAXIES 

Two kinds of local effects in galaxy clusters may produce an anisotropy and, conse¬ 
quently, a polarization in the CMB crossing these structures: the first class consists of 
gravitational effects, due to the distortion of the Hubble flow in the region of a massive 
clusterl; the latter one is due both to the inverse Compton scattering of the photons on the 
plasma electrons and to a Doppler-like shift, connected with the motion of the structure. 
In the following I will recall in some detail all these effects, pointing out the local CMB 
anisotropy they induce in order to determine the level of linear polarization generated by 
the Thomson scattering. 

a) Gravitational effects: 

1 ) Moving Gravitational lenses: 

A massive object always influences the radiation that propagates in the surrounding 
region, gravitationally distorting its path; in particular, a source observed behind a massive 
structure appears slightly different in frequency and flux if compared to a similar one 
located in a different position: that represents a typical lensing phenomenon. 

A similar, but less known, effect is produced when the object representing the grav¬ 
itational lens is in motion across the line of sight (Birkinshaw 1989); in this case, the 
lens produces a) an anisotropy in the CMB crossing the lens and b) a redshift difference 
between the multiple images of a background object seen through it. 

These effects have been discussed in (Birkinshaw 1989; Birkinshaw & Gull 1983); here 
we are interested in the production of an anisotropy in the CMB, a phenomenon that 
depends on the transverse velocity of the lens and presents a peculiar angular dependence. 

Looking at fig. 2a, we consider the CMB as seen through a convex lens (for instance a 
galaxy cluster), moving with a velocity v which forms an angle a with the line of sight: 5 
is the deflection angle of the radiation and 6 is the angle between the deflected direction of 
propagation and the line of sight (see fig. 2b) (Birkinshaw & Gull 1983). In the lens frame, 
the frequency of the radiation is unchanged: on the contrary, for a stationary observer in 
the frame of the original source, the observed frequency of the deflected light depends on 


4 



v and 5 as follows (Birkinshaw & Gull 1983): 


An 

v 


7/3 5 sin a cos </>; 


(3.1) 


here f3 and 7 respectively are v/c and the Lorentz factor, while <i> is the angle between 
the lens velocity vector and the emerging ray projected onto the sky (see fig. 2b). The 
deflected photons have different energies and different brightness in comparison with the 
undeflected ones: as a result, the corresponding change in the brightness temperature of 
the radiation field is (Birkinshaw & Gull 1983): 


AT 

~T~ 


( 0 , 0 , a) 


87 /3GM 0 sin a 
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COS 4 > ■ 
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e 2 \ 3/2 
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(3.2) 


eq. (3.2) holds for 9 < 9 0 ; M 0 is the lens mass (that, for a galaxy cluster, is about 
10 15 Mq), L is the average cluster radius (L ~ 3 Mpc ) and 9$ is given by the ratio L/d 
where d is the cluster distance. 

From eq. (3.2) one can remark that the CMB anisotropy and the induced polarization 
could be relevant for cluster of galaxies having large transverse peculiar velocities: in fact, 
in a rich cluster of galaxies, the deflection angle due to the lensing is relatively large 
(~ 1 arcmin ), while the transverse velocity may be larger than 1000 km s -1 (Dressier 
1987). 


2) Effects due to the gravitational bound: 


The decoupling of a massive, gravitationally bounded structure from the local Hubble 
flow produces a redshift in the radiation that crosses it, when compared to the one that 
passes far away; in fact, the metric of the massive object changes in a different way as 
regards the one of the surrounding Universe. As a result, a monochromatic source observed 
behind a massive structure will be observed at a shifted frequency compared to a similar 
source placed in a different position. 

This redshift effect in a variable gravitational field was theoretically studied by(Rees 
& Sciama 1968; Dyer 1976; Nottale 1984) in the framework of a “Swiss cheese” model: 
this model starts from a zero-pressure Robertson-Walker Universe where one removes a 
comoving sphere of dust and places a clump of the same mass at the centre of the “hole” 
thus obtained. In this way, the optical effects of the massive object are influenced by 
the surrounding void in the embedding Universe: that is due to the localized changes of 
the space-time metric in the region where the inhomogeneity is present and where the 
distortion of the Hubble flow is significant. 

In this model, a cluster is represented by a high density Friedmann solution of the 
Einstein equations, separated by an empty Schwarzschild zone from the surrounding, 
lower density Friedmann Universe (Nottale 1984); the radiation crossing both the empty 
Schwarzschild zone and the massive central structure suffers two effects (Dyer 1976): a) a 
redshift effect, implying that the parameters 1 + z = v/vq and 1 + £ = R(Tq)/R(T) (i? 
is the scale factor) do not coincide behind the cluster; b) a time-delay effect, due to the 
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fact that the radiation crossing the inhomogeneous structure follows a longer path. For 
the CMB photons the temperature anisotropy resulting from these effects reads (Nottale 
1984): 


AT 

~T~ 
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(3.3) 


where q 0 is the deceleration parameter, L is the cluster size, pc and po are the average 
density of the cluster and the one of the background Universe and finally II c and Hp are 
the Hubble constants for the cluster and for the background Universe at the present time. 
The ratio ~pl/Po is about equal to 10 in the case of an expanding cluster, while it is about 
1000 for a collapsing structure (Nottale 1984). 


3) Effects due to the gravitational collapse or expansion: 


A phenomenon very similar to the one above discussed is produced when the change 
of the metric is due to the collapse or to the expansion of the massive object (Nottale 
1984): a rich cluster may both contract (H c /Hp < 0) or expand ( H c /Hp > 0); however, 
the observations seem to support the possibility of a collapse (Capelato et al. 1982; de 
Vaucouleurs 1982). 

The resulting temperature anisotropy of the CMB photons is (Nottale 1984): 

a) for a contracting cluster (. H c /Hp < 0): 
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b) for a slightly expanding cluster (0 < H c /Hp < 1): 
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c) for a rapidly expanding cluster ( H c /Hp > 1): 
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(3.4) 


(3.5) 


(3.6) 


These gravitational effects seem to dominate the scattering ones; that is surely true 
for very massive clusters (M > 10 19 Mg), as proved in (Rees & Sciama 1968); in the 
following I will prove their relevance also for lower masses (M rsj 10 14 - 10 15 Mq). 


b) The Sunyaev-Zel’dovich and the Doppler effects: 


The Sunyaev-ZePdovich effect (SZ) (Sunyaev & Zel’dovich 1972) consists in a char¬ 
acteristic distortion of the CBM photons spectrum due to the inverse Compton scattering 
of the radiation on the hot electron gas (T ~ 10 8 K ) present in the clusters (Lea et al. 
1973). That produces a shift in the CMB spectrum towards smaller wavelengths near the 
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blackbody peak at A ~ 1 mm and a consequent decrement of the radiation intensity in 
that region; correspondingly, one observes an increment of the intensity at submillimetric 
wavelengths. 

This effect has been studied at centimeter wavelengths in many clusters, as Abell 2319 
(White & Silk 1980) Abell 576 (Boynton et al. 1982), (White & Silk 1980), Abell 2218 
(Boynton et al. 1982; Klein et al. 1991; Jones et al. 1993), Abell 665 and C10016+16 
(Birkinshaw et al. 1984); a significant evidence of the SZ distortion has probably been 
obtained for Abell 2218 (Boynton et al. 1982). In fact, the properties of the core gas in 
A576 and in some other clusters desagree with the theoretical expectations for the inverse 
Compton scattering process and therefore the apparent CMB decrement could be due to 
some spurious effects (White & Silk 1980). Anyway, the presumable average magnitude of 
the temperature decrement is in the range 3 x 10“ 5 -=- 10“ 3 K° (White & Silk 1980); for 
Abell 2218 Klein et al. found AT = —0.21 ± 0.04 mK° at A = 1.2 cm (Klein et al. 1991), 
while Jones et al. (Jones et al. 1993) found AT ~ 0.5 mK°. 

From a theoretical point of view, the evaluation of the magnitude of the SZ effect 
depends on the cluster model one assumes, in particular on the X ray emissivity, on the 
cluster size and mass, on the form of the gravitational potential and on the physical con¬ 
ditions characterizing the core gas. Generally, one supposes the hot gas is in an isothermal 
and quasi-hydrostatic equilibrium and the cluster is a stationary and spherically symmetric 
structure (White & Silk 1980). 

The X-ray emission mainly comes from the bremsstrahlung processes and it is given 
by (White & Silk 1980): 


A(T) = 3.2 x 10 


-24 


In (1 ± 
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Ei T E 2 
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exp(—Ei/T) — exp(—E 2 /T) 


ergs cm 3 s 1 , 


(3.7) 


where the passband is (T/, E 2 ) and T is the gas temperature, all expressed in KeV. The 
emissivity profile is given by (White & Silk 1980): 


L(r) = nl(r) A (T) 


nj(r) y 
2 r c 



(3.8) 


where Zq and r c are the central surface brightness and the X-ray core radius, known from 
the observations. 

Finally, in an isothermal model for the gas (Boynton et al. 1982; Klein et al. 1991) 
the density profile is given by the following formula: 


n e 


r 2 3 n/2 

0 + 7I> 


(3.9) 


where n ej o is the central electron density, equal to (3.4 ± 1.5) x 10 3 cm 3 , and n is the 
power law index, equal to 0.5 ±0.1 (Klein et al. 1991). 
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For a gas whose properties are given by eqs. (3.6)-(3.9), the microwave decrement at 
the centre of the cluster in the Rayleigh-Jeans limit is (White & Silk 1980): 


ATm 

Tm 


2 k(JT 
m e c 2 


(2r c l 0 ) 1/2 



T(x)(A[T(x)]) 


(1 + x 2 )- 3/i dx, 


(3.10) 


where (7t = 6.65 x 1CD 25 cm 2 is the Thomson scattering cross section and the microwaves 
temperature has been called Tm in order to distinguish it from the gas temperature T; 
finally, x = hu/kT. 

In an isothermal model with T = 12.5 KeV the integration of eq. (3.10) gives (White 


& Silk 1980): 
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Tm 


-7.4 


kar 

m e c 2 


(r c l o) 1/2 T [A(T)]- 1 / 2 ; 


(3.11) 


eq. (3.11) is the final expression I will assume for the SZ temperature decrement. 


An additional decrement is due to the peculiar motion of the cluster (Sunyaev & 
Zel’dovich 1972 and 1980): that is just a Doppler-like effect and it causes a change of the 
radiation intensity and of the CMB temperature depending on the radial component v r of 
the peculiar velocity of the cluster (Sunyaev & ZePdovich 1980): 


A Tm _ v r r 
Tm c 


(3.12) 


A I_ = x exp(x) v r r 

I ( exp(x ) - 1) c ’ 1 ’ 

here r is the optical depth, equal to r = n e (r)aT dr. 

While the thermal SZ effect gives the same perturbation both for the temperature and 
the intensity, in this case the expressions for A T/T and A I/I are different: in particular, 
the perturbed intensity only (eq. (3.13)) depends on the frequency. This effect should be 
dominant in galactic halos, where the gravitational and thermal contributions probably 
are negligible. 


4. THE CMB POLARIZATION IN A2218, A576, A2163. 

Now I can apply the formalism developed in sects. 2 and 3 to study the CMB polar¬ 
ization in some observed clusters. 

The clusters considered in this analysis were studied in many observations and their 
properties are quite well determined; moreover, they are possible candidates for which the 
observation of the Sunyaev-Zel’dovich decrement has been claimed. 

In tabs, la, lb and lc I resume the physical properties of these clusters useful for my 
calculation of the CMB polarization. 

4.1 THE CLUSTERS A2218, A576, A2163. 


Cluster A2218 



Following the original Abell catalogue, A2218 is classified as a richness class 4, distance 
class 6 cluster. Subsequently, its redshift has been fixed to 0.174 (Le Borgne et al. 1992) 
and its class fixed as Bautz-Morgan II (Leir & Van den Bergh 1978); it contains a large 
cD galaxy and many weak radio sources (Andernach et al. 1988); a review of the present 
characteristics of A2218 can be found in (Boynton et al. 1982; Klein et al. 1991). 

Cluster A576 

A576 is a moderately rich cluster, classified as Abell richness class 1, distance class 2 
and Bautz-Morgan class III (Leir & Van den Bergh 1978); it presents a central condensation 
region but none dominant galaxy; the cluster contains a high proportion of SO galaxies 
(Melnick & Sargent 1977), it is surrounded by an extensive X-ray emitting halo (Forman 
et al. 1978) but it is a quite weak X-ray source. A review of its properties can be found 
in (White & Silk 1980). 

Cluster A2163 

A2163 is a rich cluster of Rood and Sastry class I (Struble & Rood 1987) ( correspond¬ 
ing to NAbell — H9 galaxies) having a redshift z — 0.201. The spectroscopic X-ray obser¬ 
vations, performed with the GINGA satellite (Arnaud et al. 1992) showed an exceptionally 
high temperature and X-ray luminosity (kT — 13.91^'q KeV and Lx — 6.0 x 10 45 erg/s 
in the band 2 — 10 KeV, while the usual values for other clusters are respectively lower 
than 9 KeV and 2 x 10 45 erg I see). Due to the fact that the SZ effect is proportional to 
the product of the electron density and of the temperature, A2163 surely represents a very 
promising candidate for the observation of the SZ microwave decrement. A review of the 
present experimental knowledge of A2163 can be found in (Elbaz et al. 1995). 

4.2 THE CALCULATION OF THE POLARIZATION 

Before calculating the CMB polarization induced by the previously discussed anisotropies 
we should firstly fix a suitable reference frame to work. If we look at figs. 2a and 2b, a 
possible choice consists in fixing the origin of the frame in coincidence with the observer 
on the earth; the z axis is oriented towards the centre of the cluster, along the line of sight. 
For simplicity, I assume an approximate spherical symmetry for the clusters, a hypothesis 
in most cases well confirmed by the observations. 

I introduce the usual polar coordinate system (r, 6, </>) and I define 6q = L/d, where L is 
the cluster radius and d is its distance from the observer. Then, I write the transfer equation 
(2.2) in this particular reference frame. Before the scattering Thomson the anisotropic 
radiation is unpolarized: the perturbed intensity reads 

[AI]tot — [A/(0, 4>, a)}Lens + [A/(r)] Grav.b.+coll./exp. + [A I(T)\sz + [AI(u, V r )]Doppl, 

(4.2.1) 

i.e. it is given by the sum of the various contributions previously discussed. The incoming 
radiation in the transfer equation is represented by a vector ni n containing the perturbed 
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intensity: 


n m = (A//2, A//2, 0). (4.2.2) 

(as usual, I neglect the Stokes parameter V because no circular polarization is produced). 
After the Thomson scattering on the freee eleectron plasma, the anisotropic radiation 
acquires a linear polarization (Rees 1968) and the vector expressing its status is: 

n SC att = (AI h AI r , AU). (4.2.3) 


Then, the transfer equation in polar coordinates reads: 

dn sca tt 1 dn sca tt ro a , ^ 1 dn sca tt r , n 

3 — - + - • ——— 2 tana — cot 0 + - ■ ——— tan® — cot® = 

or r dO r dtp 


— (J'J 1 — 
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-27r 


d/A d(f)' P(n, n', 4>, (f)') n in ; (4.2.4) 


'MO 


(/i 0 = cos 0 o! the primed angular variables refer to the incoming radiation). 

Eq. (4.2.4) holds on fixed-time hypersurfaces: the reference time t corresponds to the 
observed redshift of the cluster. The variable 6 in the integration spans over the range 

( 0 , 0q). 


4.3 THE RESULTS 


The CMB polarization is obtained by solving eq. (4.2.4): due to the axial simmetry 
of the problem, this expression simplifies, because the term containing the 0 derivative 
vanishes. 

Substituting the vectors (4.2.2) and (4.2.3) into eq. (4.2.4), I obtain a system of three 
partial differential equations in the perturbed parameters A/;, A I r , AU, whose solution 
represents the final polarization status of the CMB in the cluster. The system reads as 
follows: 
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(4.3.1a) 

- AI r ] , 
(4.3.16) 

- AU]. 


(4.3.1c) 

In eqs. (4.3.1a), (4.3.1b) and (4.3.1c) I called R\, R 2 and R 3 the scattering integral that 
appears in eq. (4.2.4), while R(t) is the scale factor of the Universe; I separate the different 
contributions to the anisotropy in the following way: 



tot 


= F^O', (j)', a) + F 2 (r) + F 3 (u,T), 


(4.3.2) 
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where F\ is given by eq. (3.2) and represents the lensing contribution (the primed angular 
variables must be used because we are referring to the incoming radiation) F 2 is given 
by eq. (3.3) plus one of the eqs. (3.4), (3.5) or (3.6), depending on the cluster evolution 
(gravitational bound + collapse or expansion term) and finally F 3 is given by the sum of 
eq. (3.11) (thermal SZ effect) and eq. (3.13) (Doppler effect); then the scattering integrals 
read (j — 1, 2, 3): 


r>2i r 




/ / dfx d(j)' {a jl +a j2 ) F 1 {9 ', </>', a) + F 2 (r) + F 3 (v,T) 

J 0 J cos do L 

(4.3.3) 

In eq. (4.3.3), a 3 \ and cij 2 are the matrix elements of the Chandrasekhar scattering matrix, 
whose form can be found in Appendix. 

A possible way to perform the calculation consists in searching a recursive solution of 
eqs. (4.3.1a), (4.3.1b) and (4.3.1c): I expand the perturbations in Legendre polynomials 
as follows: 

All = — ^2 ai p i( cos ®), 


AL = 


^2 bi p i ( cos6 )» 


AU = — V' Pz(cosd), 

nr* ( ^ 


(4.3.4a) 

(4.3.46) 

(4.3.4c) 


where Iq is the unperturbed intensity, given by: 


In = 


exp(hv/ksT) — 1 ’ 


(4.3.5) 


here hV ~ 10 -13 GeV and ksT ~ 2.35 x 10 -13 GeV. Then, it is possible to calculate the 
coefficients of the expansion simply by solving the differential equations (4.3.1a), (4.3.1b) 
and (4.3.1c); the recursive calculation can be stopped at l = 2, because the higher order 
terms in the polynomials are negligible. Finally, the polarization degree P is 




All ~ AI r + AU 
(Ah + AI r ) 


(4.3.6) 


In tabs. 2a, 2b and 2c I show a comparison of the different contributions F±, F 2 and F : > to 
the local CMB anisotropies for the considered clusters and in the cases of a contracting, 
slightly expanding or totally expanding cluster; the data listed in these tables have been 
calculated at some reference values for the polar coordinates and for an angle a = 7t/2 in 
order to maximize the lensing contribution. In tab. 3 I list the values of the coefficients 
which appear in the expansions (4.3.4a), (4.3.4b) and (4.3.4c) of the perturbed Stokes 
parameters for the three clusters: I show the values obtained for fixed (r, 6 , </>) in order 
to give an idea of the size of these coefficients. Finally, in figs. 3, 4 and 5 I plotted the 
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average polarization degree in the clusters A2218, A576 and A2163 as resulting from this 
calculation. 

5. DISCUSSION AND CONCLUSIONS 

Looking at eqs. (4.3.1a), (4.3.1b) and (4.3.1c), we remark that the different evolution 
of the perturbed Stokes parameters originates in the scattering integrals, as an effect of 
the dissimilarity of the functions R,\ , R -2 and R 3 . 

The angular dependence of the matrix elements cij 1, aj2 (see eqs. (A.l), (A. 2 ) and 
(A.3) of the Appendix) and of the gravitational lensing term iq plays a fundamental 
role in determining the Stokes parameters evolution and, therefore, in producing a CMB 
polarization; for instance, the R 3 scattering integral vanishes when integrated on 0 f : as 
a result, all the coefficients <q of the expansion (4.3.4c) are equal to zero and the overall 
contribution of A U to the polarization cancels. A linear polarization results from the 
scattering contributions Ri and /A: the difference existing between these functions causes 
an unequal evolution of the perturbed left and right intensities and therefore an unvanishing 
term Al) — A/ r | in eq. (4.3.6). 

In tabs. 2a, 2b and 2c, I considered both the possibilities of a collapse or an expan¬ 
sion for the clusters: we can remark that the most important contribution to the CMB 
anisotropy and polarization is due to the gravitational collapse of the clusters (7q> term); 
the high densities attained in this case strongly enhance the importance of the gravitational 
effects but, anyway, also for an expanding cluster F2 dominates. Moreover, the SZ effect 
(iq contribution) dominates over the lensing one (T\ contribution), while the Doppler shift 
is totally negligible at the considered frequencies; note however that the lensing contribu¬ 
tion slightly influences the angular behaviour of the scattering integral, by introducing an 
additional cos </>' dependence. 

In tab. 3 I listed the coefficients of the expansion in Legendre polynomials, obtained by 
solving the differential equations (4.3.1a), (4.3.1b) and (4.3.1c): the recursive calculation 
can be stopped at the order l — 2 because the coefficients vanish very fast. 

These results confirm that the gravitational effects are important in determining the 
characteristics of the CMB photons crossing the galaxy clusters, particularly in the case 
of very massive and dense structures. In figs. 3, 4 and 5 I showed the behaviour of the 
average polarization degree in galaxy clusters as a function of the polar angle 9: due to the 
smallness of 9 max — 9 q — L/d, this dependence is very weak, indicating that an observation 
performed towards the clusters centres should be able to measure the CMB polarization 
induced by these effects. Possible fluctuations in the signal might be attributable to the 
presence of some central radiosources which introduce a serious source of noise. 

The average polarization degree for A2218 and A2163 is comparable, while for A576 
it is smaller by a factor 10: that is probably due to the z dependence of the gravitational 
bound effect. 

No relevant differences are found between the polarization in expanding or in collapsing 
clusters: that is due to the fact that a change in the function F2 affects in the same way 
the values of the left and right perturbed intensities, thus leaving unchanged the quantity 
A// - A/ r . 
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Summarizing, in this paper I studied the CMB polarization in clusters of galaxies 
induced by local anisotropies: by expanding the Stokes parameters in Legendre polynomials 
and by solving the Boltzmann transfer equation through a recursive method, I proved that 
the gravitational effects have some importance: they indeed influence the properties of the 
CMB photons crossing the galaxy clusters. The cancellation of some terms in the Thomson 
scattering integral makes A/; — AI r \ different from zero, while the resulting perturbation 
in the Stokes parameter U vanishes. 

In such a way, the Thomson scattering of the locally anisotropic CMB radiation (the 
anisotropy being due both to the scattering and to the gravitational effects) produces 
an appreciable polarization in clusters having a mass M > 10 14 -j- 10 15 Mq An accurate 
observation towards the clusters centres might, in principle, reveal such a phenomenon, 
provided that one is able to distinguish the contribution due to the possible presence of 
strong radiosources. 
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APPENDIX 


The Chandrasekhar scattering matrix (Chandrasekhar 1960) is a 4 x 4 matrix having 
the following structure: 


/ On 

012 

O13 

0 

\ 

021 

022 

023 

0 


031 

032 

033 

0 


\ 0 

0 

0 

CL 44 

/ 


The presence of only one matrix element different from zero in the position (4,4) assures 
the Stokes parameter V is totally decoupled from the other ones: as a result, the Thomson 
scattering of the anisotropic radiation does not produce a circular polarization. 

When we apply the scattering matrix to the vector (4.2.2), representing the incoming 
radiation, we need to know the following linear combinations of the matrix elements a^: 
for eq. (4.3.1a): 


an + ai2 


3/j 2 /j 2 — /r 2 — 2/i 2 + 2 + 4 y/l — n 2 \/l — /r' 2 mi cos(4>' — </>) + 


for eq. (4.3.1b): 


+ /j 2 /j 2 cos 2 (cf)' — (f)) — ft 2 cos 2 ((f)' — (f>) ; 


(A.l) 


0,21 + a 22 


1 + fi 2 + cos 2 (<f>'— <f>) (1 — n 2 ) ; 


(T.2) 


for eq. (4.3.1c): 


3 

«31 + O32 — - 


-2/rVl-M 2 V 7 ! — l 1 ' 2 sin {(f)' — (f)) + sin 2 ((f) 1 — (f)) (// — /j/i 2 ) 


(A. 3) 


The (f)' dependence of eqs. (A.l), (A.2) and (A.3) produces the cancellation of some terms 
in the scattering integrals; the anisotropy due to the lensing effect (see eq. (3.2)) also 
contains a factor cos </>', thus, combining these dependences, I obtain from eq. (4.3.3) 
Rs — 0 and Ri R 2 . This difference between Ri and R 2 is finally responsible for the 
creation of a local polarization in clusters. 
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Tab. la: Properties of the cluster Abell 2218 (Boynton et al. 1982; Klein 
et al. 1991). 


CLUSTER PROPERTIES 

DATA 

z 

0.174 

d 

1060 ± 810 Mpc 

L 

~ 3 Mpc 

Oo 

2.83 x 10“ 3 

r c 

0.22 ± 0.06 Mpc 

lo 

1.60 x 10 -4 erg cm~ 2 sec~ 4 

n 0 

(4.2 ± 1.8) x 10“ 3 cm~ 3 

X 

0.438 {y = 25 GHz) 

Mt 

(7.8 ± 1.4) x 10 14 Mq 

P 

0.02 

o 

1400 ± 200 Km/s 

n 

0.5 


Note: d = cluster distance; L = cluster radius; 9 0 = L/d ; r c = X-ray 
core radius; Iq = central surface brightness; no = central electron density; 
x — hv/kTr ; Mt = total cluster mass; (5 = v c i u /c ; a = line of sight velocity 
dispersion of the galaxies, n = power law index of the electron density profile. 
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Tab. lb: Properties of the cluster Abell A576 (White & Silk 1980) 


CLUSTER PROPERTIES 

DATA 

z 

0.039 

d 

236 Mpc 

L 

~ 3 Mpc 

do 

0.013 

r c 

0.24 Mpc 

lo 

3.8 x 10~ ;> erg cm~ 2 sec~ x 

n 0 

4.6 x 10~ 3 cm -3 

X 

0.18 {y = 11 GHz) 

Mt 

(2.9 ± 0.5) x 10 15 M q 

P 

0.02 

a 

1124 Km/s 

n 

0.64 
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Tab. lc: Properties of the cluster Abell 2163 (Struble & Rood 1987; 
Arnaud et al. 1992; Elbaz et al. 1995) 


CLUSTER PROPERTIES 

DATA 

z 

0.201 

cl 

1267 Mpc 

L 

~ 3 Mpc 

Oo 

2.37 x 10“ 3 

r c 

0.305 ± 0.019 Mpc 

lo 

2.16 x 10 -3 erg cm ~ 2 sec “ 1 

n 0 

6.8 x 10 -3 cm~ 3 

X 

0.18 {v = 11 GHz) 

Mt 

(4.6 ± 0.4) x 10 15 M q 

P 

0.02 

o 

U nknown 

n 

0 62 +0 ' 05 
u - dz -0.02 
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Tab. 2a: Different contributions to the CMB anisotropy and polarization 
in A2218. 


Ei 

E 2 (GB) 

E 2 (CO/EXP) 

e 3 


28.690 

-2484.92 (CO) 


1.63 x 10“ 6 

0.116 

-0.4720 (SE) -3.923 x 10“ 4 



0.116 

0.142 (EXP) 



F \: Lensing contribution; 

F 2 : Gravitational bound contribution (GB); Gravitational collapse (CO) 
or expansion (Slight expansion (SE), Expansion (EXP)) contributions; 

F 3 : Sunyaev-Zel’dovich contribution (SZ); the Doppler shift contribution 
is negligible. 

The listed values correspond to the following choice of the cluster coordi¬ 
nates and properties: r = 1063 Mpc, 0 = 2.83x10 3 rad , <f> = 0.785 rad , a = 1.570 rad. 
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Tab. 2b: Different contributions to the CMB anisotropy and polarization 
in A576. 


Ti 

F 2 (GB) 

F 2 (CO/EXP) 

f 3 


0.299 

-125.62 (CO) 


5.89 x 10“ 6 

1.047 x 10“ 3 

-5.361 x 10“ 3 (SE) 

-3.183 x 10“ 5 


1.047 x 10“ 3 

1.342 x 10“ 3 (EXP) 



F \: Lensing contribution; 

F 2 : Gravitational bound contribution (GB); Gravitational collapse (CO) 
or expansion (Slight expansion (SE), Expansion (EXP)) contributions; 

F 3 : Sunyaev-Zel’dovich contribution (SZ); the Doppler shift contribution 
is negligible. 

The listed values correspond to the following choice of the cluster coordi¬ 
nates and properties: r = 239 Mpc, 9 = 0.013 rad, (f> — 0.785 rad, a — 1.570 rad. 
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Tab. 2c: Different contributions to the CMB anisotropy and polarization 
in A2163. 


Ei 

E 2 (GB) 

F 2 (CO/EXP) 

e 3 


49.773 

-3546.93 (CO) 


9.56 x 10“ 6 

0.206 

-0.804 (SE) -5.411 x 10“ 5 



0.206 

0.250 (EXP) 



F \: Lensing contribution; 

F 2 : Gravitational bound contribution (GB); Gravitational collapse (CO) 
or expansion (Slight expansion (SE), Expansion (EXP)) contributions; 

F 3 : Sunyaev-Zel’dovich contribution (SZ); the Doppler shift contribution 
is negligible. 

The listed values correspond to the following choice of the cluster coordi¬ 
nates and properties: r = 1270Mpc, 6 — 2.37xl0~ 3 rad , </> = 0.785 rad , a — 1.570 rad. 
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Tab. 3: The coefficients of the expansion in Legendre Polynomial of the 
perturbed Stokes parameters: a; refers to A/;, bi to A I r while q, referring to 
A U, are zero for any /; the values have been calculated for a fixed r, respectively 
r=1063 Mpc for A2218, r=239 Mpc for A576 and r=1270 Mpc for A2163; I 
chose also 6 = 9q and — 7r/4. 


CLUSTER 

ai 

bi 


a 0 = -1.206 x 10+ 25 

b 0 = -7.654 x 10+ 24 

A2218 

a x = -7.638 x 10+ 17 

bx = -7.636 x 10+ 17 


tt2 ~ 0.0 

62 ~ 0.0 


a 0 = -2.801 x 10+ 24 

b 0 = -2.750 x 10+ 24 

A576 

ax = -1.717 x 10+ 17 

bx = -1.738 x 10+ 17 


Cl2 ~ 0.0 

&2 ~ 0.0 


a 0 = -1.440 x 10+ 25 

b 0 = -6.895 x 10+ 24 

A2163 

ax = 9.126 x 10+ 17 

bx = 4.563 x 10+ 17 


tt2 ~ 0.0 

62 ~ 0.0 
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FIGURE CAPTIONS 


Fig. 1: 

Fig. 2a: 

Fig. 2b: 

Fig. 3 
Fig. 4 
Fig. 5 


The polarization ellipse (see also (Chandrasekhar 1960), p.26); the princi¬ 
pal axes of the ellipse form the angles y, X + 7r /2 with the direction l] ip is 
the angle formed by the generical vibration direction with /. The Stokes 
parameters are defined as a function of 1), / r , y and f3 (the tangent of (3 is 
the ratio of the axes of the ellipse traced by the end point of the electric 
vector) as in (Chandrasekhar 1960). 

A moving lens produces a change in the brightness of an isotropic radiation 
field, proportional to the transverse velocity of the lens: the deflection 
causes a slight decrease in the photon energy if 6 > 0 and a slight increase 
if 6 < 0; in the Rayleigh-Jeans part of the CMB spectrum this effect appears 
as a brightness increase or decrease. The figure is taken from (Birkinshaw 
& Gull 1983). 

The configuration of the problem in the frame of the observer; the lens has 
a velocity v in the (x, z) plane, while the angles (6, <p) describe the direction 
of the deflected photon and respectively correspond to the angle between k 
and the 2 axis and to the angle between the projected k and the projected 
v. The figure is taken from (Birkinshaw & Gull 1983). 

The average polarization degree as a function of the polar angle 6 in the 
cluster A2218. 

The average polarization degree as a function of the polar angle 6 in the 
cluster A576. 

The average polarization degree as a function of the polar angle 6 in the 
cluster A2163. 
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